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ON INCOMPRESSIBLE NAVIER-STOKES DYNAMICS: A NEW 
APPROACH FOR ANALYSIS AND COMPUTATION 

JIAN-GUO LIU, JIE LIU, AND ROBERT L. PEGO 


Abstract. We show that in bounded domains with no-slip boundary condi¬ 
tions, the Navier-Stokes pressure can be determined in a such way that it is 
strictly dominated by viscosity. As a consequence, in a general domain we can 
treat the Navier-Stokes equations as a perturbed vector diffusion equation, 
instead of as a perturbed Stokes system. To illustrate the advantages of this 
view, we provide a simple proof of the unconditional stability of a difference 
scheme that is implicit only in viscosity and explicit in both pressure and con¬ 
vection terms, requiring no solution of stationary Stokes systems or inf-sup 
conditions. 


1. Introduction 

Consider the Navier-Stokes equations for incompressible fluid flow in f 1 with 
no-slip boundary conditions on T := dCl: 

(1) d t u + u-Vu + Vp = vAu + f in Q, 

(2) V • u = 0 in Q. 

(3) u = 0 on r. 

Here u is the fluid velocity, p the pressure, and v is the kinematic viscosity coeffi¬ 
cient, assumed to be a fixed positive constant. Pressure plays a role like a Lagrange 
multiplier to enforce the incompressibility constraint, and this has been a main 
source of difficulties. In this paper we will describe some of the main results of 
which indicate that the pressure can be obtained in a way that leads to 
considerable simplifications in both computation and analysis. 

A standard way to determine p is via the Helmholtz-Hodge decomposition. We 
can rewrite m as 

(4) d t u + V(u -Vu — / — pAu) = 0. 

where V is the standard Helmholtz projection operator onto divergence-free fields. 

In this formulation, solutions formally satisfy c^(V • u) = 0. The dissipation in 
0 appears degenerate due to the fact that V annihilates gradients, so the analysis 
of ® is usually restricted to spaces of divergence-free fields. But alternatives are 
possible if the pressure is determined differently when the velocity field has non-zero 
divergence. Instead of 0 , we will consider the unconstrained formulation 

(5) d t u + V(u-Vu- f-vAu) =zA7(V-u). 

There is no difference as long as V • u = 0, of course. But we argue that 0 enjoys 
superior stability properties, for two reasons. The first is heuristic. The incom¬ 
pressibility constraint is enforced in a more robust way, because the divergence of 
velocity satisfies a diffusion equation with no-flux boundary conditions. Naturally, 
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if V • u = 0 initially, this remains true for all later time, and one has a solution of 
the standard Navier-Stokes equations JU ©. 

The second reason is more profound. To explain, we recast © in the form © 
while explicitly identifying the separate contributions to the pressure term made 
by the convection and viscosity terms. Using the Helmholtz projection operator V, 
we introduce the Euler pressure p E and Stokes pressure ps via the relations 

(6) V(u-Vu- f) = uA7u- f + \7p E , 

(7) V{—Au) = — Au + V(V • u) + Vp s - 
This puts © into the form © with p = p E + vp s : 

(8) dtu + u -Vu + Vp E + uVp s = vAu + f. 

Since the Helmholtz projection is L 2 -orthogonal, naturally the Stokes pressure from 
© satisfies 

(9) / |Vp s | 2 < / |Aii| 2 if V • u = 0. 

Jn Jn 

In fact, it turns out that the Stokes pressure term is strictly dominated by the 
viscosity term, regardless of the divergence constraint. The estimate contained 
in the following theorem from elei is key to a new treatment of Navier-Stokes 
dynamics in bounded domains. 

Theorem 1. Let H C (N > 2) be a connected bounded domain with C 3 
boundary. Then for any e > 0, there exists C > 0 such that for all vector fields 
u £ H 2 fl Hq(Q,S. n ), the Stokes pressure ps determined by © satisfies 

(10) [ \V Ps \ 2 < P f |Am| 2 + C [ |Vm| 2 . 

JQ J £1 J £1 

where /3 = | + e. 

In this paper we will prove a slightly weaker version of this result that has a 
simpler proof. Namely, we will show that the estimate ©1 holds for some ft < 1. 
The full result that fill holds with any ft > 2/3, independent of the domain or 
the space dimension, is proved in EES using differential geometry to establish 
sharp integrated Neumann-to-Dirichlet estimates for functions harmonic near the 
boundary of fh 

Due to Theorem© we can treat the Navier-Stokes equations in bounded domains 
simply as a perturbation of the vector diffusion equation dtu = vAu, regarding both 
the pressure and convection terms as dominated by the viscosity term. To begin 
to see why, recall that the Laplace operator A: H 2 (Ll) fl -ffg(H) —> L' 2 (fl) is an 
isomorphism, and note that Vps is determined by Au via 

(11) Vp s = (7 - V - Q) Au, Q := VV ■ A" 1 . 

Equation © can then be written 

dtu + V(u-Vu — /) = v{V + Q)Au 

(12) = vAu — 1/(1 — V — Q) Au. 

Theorem ©will allow us to regard the last term as a controlled perturbation. 

This approach should be contrasted with the usual one that regards the Navier- 
Stokes equations as a perturbation of the Stokes system. To show the advantage 
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of this point of view, we will sketch the proof from CLP 1 of unconditional stability 
of a simple time-discretization scheme with explicit time-stepping for the pressure 
and nonlinear convection terms and that is implicit only in the viscosity term. 

The discretization that we use is related to a class of extremely efficient nu¬ 
merical methods for incompressible flow ITI1 IPel I.TLI IGuSI . Thanks to the explicit 
treatment of the convection and pressure terms, the computation of the momentum 
equation is completely decoupled from the computation of the kinematic pressure 
Poisson equation used to enforce incompressibility. No stationary Stokes solver is 
necessary to handle implicitly differenced pressure terms. For three-dimensional 
flow in a general domain, the computation of incompressible Navier-Stokes dynam¬ 
ics is basically reduced to solving a heat equation and a Poisson equation at each 
time step. This class of methods is very flexible and can be used with all kinds of 
spatial discretization methods [.I 1 .1 . including finite difference, spectral, and finite 
element methods. The stability properties we establish here should be helpful in 
analyzing these methods. 

Indeed, we will show here that our stability analysis easily adapts to proving 
unconditional stability for corresponding fully discrete finite-element methods with 
C 1 elements for velocity and C° elements for pressure. (For additional details and 
convergence results see CEP].) It is important to note that we impose no inf-sup 
compatibility condition between the finite-element spaces for velocity and pressure. 

Below, we will also describe how the unconstrained formulation above can be 
extended to handle nonhomogeneous boundary conditions. In 11.1. FI this formula¬ 
tion was exploited to establish a new result regarding existence and uniqueness for 
strong solutions with minimal compatibility conditions for total flux through the 
boundary. For further results and more complete references and proofs, we refer to 

CLP]. 


2. Properties of the Stokes pressure 

Let G C K w be a bounded domain with C 3 boundary F. For any x G we let 
<!>(:?) = dist(a;,r) denote the distance from x to T. For any s > 0 we denote the set 
of points in Q within distance s from T by 

(13) = {x £ Q | $(a?) < s}, 

and set and r s = {f € U $(£) = s}. Since F is C 3 and compact, 

there exists so > 0 such that $ is C 3 in Q So and its gradient is a unit vector, with 
|V<&(:r)| = 1 for every x £ G So . We let 

(14) n(x) = — V<S>(x), 

then n(x) is the outward unit normal to T s = dflg for s = ^(x), and n £ 

c 2 (n S0 ,R N ). 

We let {f,g) a = fnfg denote the L 2 inner product of functions / and g in 
f2, and let || • ||n denote the corresponding norm in L 2 (fl). We drop the subscript 
on the inner product and norm when the domain of integration is understood in 
context. 

2.1. Estimates related to the Neumann-to-Dirichlet map. Taking the di¬ 
vergence and normal component of both side of 0 , we see that Stokes pressure 
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is harmonic in f l, and is determined as the zero-mean solution of the Neumann 
boundary-value problem 

(15) Ap s =0 in Ll, n ■ Vps = n • (A — VV-)S on F. 

A useful ingredient in our proof of is an integrated version of a standard 
estimate that controls tangential gradients at the boundary (these are determined 
from Dirichlet data) in terms of the normal gradient (Neumann data). 

Lemma 1. Let LI C be a bounded domain with C 2 boundary. Then there exist 
positive constants ro, C\ and Co such that for any p that satisfies A p = 0 in LI, for 
any s G (0, ro] we have 

(16) / \Vp\>< Cl J r |n.Vp| 2 , 
and furthermore, 

(17) f \V P \ 2 <Ci [ \n • Vp| 2 . 

J J Q s 

Proof: If LI has smooth boundary and s G (0, .So] is fixed, with s 0 taken as above, 
the estimate (Et is a consequence of classical elliptic theory as developed in the 
book of Lions and Magenes iLMl . applied to the Neumann problem 

(18) Ap = 0 in fig, n-Vp = g on T s . 

Under the solvability condition f r g = 0, this theory yields a map g h> p h> p|r„ 
from L 2 (r s ) —> H 3 / 2 (Ll s ) —> H 1 (T S ) that is bounded. This bounds the tangential 
part of Vp in L 2 (r s ), which is sufficient to establish (1161) . We need to verify 
that these bounds are uniform in s for domains with C 2 boundary. We defer 
details to an appendix, where we sketch a streamlined version of the arguments 
of Lions and Magenes |LM| that yield the desired bounds, based on using trace 
theorems and the Lax-Milgram lemma to solve (II81) with g G iF _1 / 2 (F s ), and 
standard elliptic regularity theory to treat g G H 1 ^ 2 (T S ), then interpolating to 
handle g G L 2 (r s ). The key to obtaining uniform estimates is to work with a fixed 
system of local boundary-flattening maps that work simultaneously for all r s - 
these are constructed using the distance function <h(a:). 

Once (ITTfl) is established, one obtains m by simply integrating m in s. 

2.2. Identities at the boundary. A key part of the proof of Theorem^ involves 
boundary values of two quantities that involve the decomposition of u = {I — 
nn t )u + nnfu into parts parallel and normal to the boundary, for which we have 
the following Lemma. The proof involves only straightforward computations and a 
density argument; please see Clp 3 for the details. 

Lemma 2. Let U C be a bounded domain with boundary T of class C 3 . Then 
for any u G H 2 (Ll,M. N ) with w|r = 0, the following is valid on T: 

(*) V • ((/ - nn t )u) = 0 in H x ' 2 { r). 

(ii) n ■ (A — VV-)(fm*it) = 0 in H~ 1 / 2 {T). 
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2.3. Identities for the Stokes pressure. Given u £ H 2 (Ll, M. N ) D 77g (SI,1^), 
using the fact that "PV(V ■ u) = 0. from 0 we get 

(19) Vp s = (7 — V)Au — VV • u = (I -V ){A - VV-)u 

Note that from the trace theorem for 77( div; S7) I GrRl l Theorem 2.5], whenever a £ 
L 2 (fl,R N ) and V • a = 0 and n • a|r = 0, then (I — V)a = 0. Thus, the Stokes 
pressure is not affected by any part of the velocity field that contributes nothing to 
ft ■ a|r where a = (A — VV-)u. Indeed, this means that the Stokes pressure is not 
affected by the part of the velocity field in the interior of f 1 away from the boundary, 
and it is not affected by the normal component of velocity near the boundary, since 
n ■ (A — VV-)(nn t u)|r = 0 by Lemma[2j 

This motivates us to focus on the part of velocity near and parallel to the 
boundary. We make the following decomposition. Fix s > 0 so 2s < min(so,r'o) 
and fix any smooth non-negative cutoff function £ : fl —► [0,1] with ( = 1 on 
Q s /2 = {.t £ | dist(ai,r) < ^s} and ( = 0 on f2\I2 s . Then we can write 

(20) u = u _l + it|| 
where 

(21) u l = + (1 — £)m, {t|| = £(I — nn t )u. 

Now u± = ( nn t )u in Q s , and with a± = (A — VV-)ttj_ we have 

(22) aj_ € L 2 (fl,R Ar ), V • a± = 0, and n-a_i_|r = 0 
by LemnraEIii). Hence (a±, V^>) = 0 for all <j> £ U 1 (H), i.e., 

(23) (7 — V)(A — VV-)it_L = 0. 

Combining this with and m proves part (i) of the following. 

Lemma 3. Let H C R w be a bounded domain with C 3 boundary, and let u £ 
77 2 (Sl,M Ar ) Cl. Letps and u\\ be defined as in \lfA) and m respectively. 
Then 

( i ) The Stokes pressure is determined by ft|| according to the formula 

(24) S7p s ={I-V)( A-VV-)t?||. 

(**) For any q £ 77 1 (H) that satisfies A q = 0 in the sense of distributions, 

(25) (Au||-Vp s ,Vg) = 0. 

(Hi) In paHicular we can let q= ps in (ii), so (Au|| — Vps, Vps) = 0 and 

(26) II Az?n || 2 = || AT?,, — Vp s || 2 + ||Vp s || 2 . 

Proof: We already proved (i). For (ii), note by Lemma 0i) we have 

(27) V • u,|| r = 0, 

so V • u\\ £ Hq (H). Using part (i), whenever q £ 77 1 (fl) and A q = 0 we get 

(28) (Vp s , Vg) = (Am|| - VV • i7||, Vg> = <Au||, Vq). 

This proves (ii), and then (iii) follows by the L 2 orthogonality. □ 
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2.4. Proof of (II OH for some (3 < 1. Let u £ R w ) n Hq (0, R n ) and define 

the Stokes pressure Vps by 0 and the decomposition u = u± + u || as in the 
previous subsection. Then by part (iii) of Lemma El we have 

(29) ||Au|| 2 = ||Au_l|| 2 +2(Am_l, Au||) + ||Ait|| - Vp s || 2 + ||Vp s || 2 - 

We will establish the Theorem with the help of two further estimates. 

Claim 1: For any e > 0, there exists a constant C e > 0 independent of u such that 


(30) (Auj_, Aft||) > -£\\Au\\ 2 -^C e \\Vu\\ 2 . 

Claim 2: There exist positive constants (3\ and C 2 independent of u such that 

(31) ||Au,| - Vp s || 2 > /3i||Vp s || 2 - C 2 ||W|| 2 . 

Proof of claim 1: From the definitions in GS> we have 

(32) A«j_ = ^nft* Au + (1 — £) Au + Ri, Au|| = £(/ — nn t )Au + f? 2 , 
where ||f?i|| + ||i? 2 || < C||Vu|| with C independent of u. Since I —fin* = (I —tin*) 2 , 

{^rtffAu + (1 — f)Au) • (£(I — nft t )Au) = 0 + £(1 — £)|(i — nn t )Au '\ 2 > 0. 

This means the leading term of (Awj_, A£t||) is non-negative. Using the inequality 
|(a,6)| < (e/C)||a|| 2 + (4C'/£)||6|| 2 and the bounds on Ri and R 2 to estimate the 
remaining terms, it is easy to obtain E3- 
Proof of claim 2: Recall that u\\ is supported in Q S: with 

(33) Att|| = £(/ — nn*) Au + R 3 
where 1111 < CHVuH. Since ft ■ (J — ftff?)Au = 0 we find 

(34) ||fJ-At?|||| n . <C 2 ||Vu|| 

with C 2 > 0 independent of u. Using ||a + 6|| 2 + ||a|| 2 > |||6|| 2 and LemmaE we 
get 

||Au|| - Vps|| 2 > [ |Vp s | 2 + f \n- (Auy - Vp s )| 2 

jqc j n a 

> [ |Vp s | 2 + ^/ \n-Vp s \ 2 -f |n-Au||| 2 
Jfil z Jn s J o s 

(35) >/3i||Vp s || 2 -C 2 ||Vu|| 2 . 

with Pi = min(l, l/(2Ci)). This establishes Claim 2. 

Now we conclude the proof of the theorem. Combining the claims with 03 , we 
get 

(36) (1 + 2e)||Au|| 2 > (1 + /3i)||Vp s || 2 - (C 2 + 2C e )||Vu|| 2 . 

Taking £ > 0 so that 2 £ < Pi yields itTTI with P = (1 + 2e)/(l + Pi) < 1. □ 
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3. Unconditional stability of time discretization with pressure 

explicit 

In this section we exploit Theorem to establish the unconditional stability of 
a simple time discretization scheme for the initial-boundary-value problem for ©, 
our unconstrained formulation of the Navier-Stokes equations. We focus here on 
the case of two and three dimensions. 

Let U be a bounded domain in R^ (TV = 2 or 3) with boundary T of class C 3 . 
We consider the initial-boundary-value problem 


(f > 0, x G fi), 
(£ > o, x g r), 
(t = o, x g Si). 


(37) dtu + u■ Vw + Vp E + pVps = vA u + f 

(38) u = 0 

(39) U = Min 

We assume {?;„ G Hq(SI,'R n ) and /G L 2 (0, T; L 2 (S2, R N )) for some given T > 0. 

Theorem Q tells us that the Stokes pressure can be strictly controlled by the 
viscosity term. This allows us to treat the pressure term explicitly, so that the 
update of pressure is decoupled from that of velocity. This can make corresponding 
fully discrete numerical schemes very efficient (see cm, also [HI El ESS])- Here, 
through Theorem [Q we will prove that the following spatially continuous time 
discretization scheme has surprisingly good stability properties: 


(40) 

(41) 

(42) 

(43) 
We set 

(44) 


u n+1 - m” 
At 


- vAu a+1 = f n -if 1 • Vtr 1 - Vp" - I/Vp£, 

Vp£ = (/-p)cr-ir-VM"), 
Vpg = (7 - V)AiT - V(V • m"), 
m"L = 0. 


rn _ 

1 ~ At 


r*(n+l) At 


iAt 


f(t)dt , 


and take u° G H 2 nTf^HjR^) to approximate Mi n in 77 q (fl, R-^). It is evident 
that for all n = 0,1, 2, ..., given ? G H 2 n Hq one can determine Vp E G L 2 and 
Vpg G T 2 from 1411) and 1421) and advance to time step n + 1 by solving (1401) as an 
elliptic boundary-value problem with Dirchlet boundary values to obtain u n+1 . 

Let us begin making estimates — our main result is stated as Theorem [2 below. 
Dot (1401) with -Au" +1 and use (1411) and ||7 — V\\ < 1 to obtain 

m" +1 || 2 


^(||ViT +1 || 2 - IIVm 71 !! 2 + ||Vm" + 1 - Vu n || 2 ) + v\\A 
<||AiT +1 ||(2||r- 


vtr 


W\\) 


(45) <^|| A? T+i||2 + l 

Z £l 


nr -tt 


Vm 71 




,'7«+ 1 ll 2 . 


IIVp^ll 2 ) 


r*n-\- 11| 2 


for any eq > 0. This gives 

■^(l|VM n+1 || 2 - ||VM n || 2 ) + (v — £ 1 )||Am 
( 46) < - (||/1 2 + ||tT • VmI 2 ) + HlVp^l 2 


£l 
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By Theorem Q for some f3 < 1 one has 

(47) HlVpsf < v(3\\AiT \\ 2 + uC/ 3 ||Vu n || 2 . 

Using this in m, one obtains 

^(||ViT +1 || 2 -||Vir|| 2 ) +(u-e 1 )(||AiT +1 || 2 - ||AiT|| 2 ) 

+ (v-e 1 -vf3)\\Au n \\ 2 

(48) < -(liril 2 + K • ViT|| 2 ) + uC^HVu n || 2 . 

£l 

At this point the pressure has been dealt with. Then, it is rather straightforward 
and standard to derive from Ladyzhenskaya’s inequalities (see [LLP] ! that 

(49) IK • VKI 2 < £2||A?T|| 2 + C||V?2 ri || 6 , 

for any £2 > 0. Plug this into © and take £ 1 , £2 > 0 satisfying v — £1 >0 and 
£ := v — £1 — v(3 — 8£2 /ei > 0. We get 

^(l|Vu” +1 || 2 - l|Vu"|| 2 ) + ( Z z-£ 1 )(||Au” + 1 || 2 - ||A?T|| 2 ) +£||Au"|| 2 

(50) < —||/"|| 2 + C||Vit"|| 6 + KsIlVu" II 2 . 

£l 

A Gronwall-type argument now leads to a simplification of the stability result in 
jLLP]: 

Theorem 2. Let LI be a bounded domain in M. N (N = 2 or 3) with C 3 boundary, 
and assume f £ L 2 (0, T; L 2 (Ll, R^)) for some given T > 0 and uP £ Hq(LI,WL n ) n 
H 2 (Cl,M. N ). Consider the time-discrete scheme 14 (V - 144V - Then there exist positive 
constants T* and C 3 , such that whenever nAt <T*, we have 

n 

(51) sup ||Vu fc || 2 + yd ||Au fc || 2 Af < C 3 . 

0<k<n fc=Q 

The constants T* and C 3 depend only upon LI, v and 

M 0 := ||Vu°|| 2 + uAt||Au°|| 2 + T \\f\\ 2 . 

Jo 

Proof: Put 

(52) *„ = ||ViZ n || 2 + ( I /-£ 1 )At||Atr , || 2 , w n = £ IIAiTH 2 , b n = ||/"|| 2 , 

and note that from iTHI) we have that Y^k =0 At < / Q T \f{t )\ 2 dt as long as nAt < T. 
Then by (l5()ll . 

(53) z n + 1 + w n At < z n + C At (b n + z n + z 2 ), 

where we have replaced max{8/£i, C , nC/ 3 } by C. Summing from 0 to n — 1 yields 

n— 1 n— 1 

(54) z n + y^uifeAt < CM 0 + CAt^2(z k + zl) =: y n . 

k—0 k =0 

The quantities y n so defined increase with n and satisfy 

(55) y n + 1 -y n = CAt{z n + z 3 ) < CAt{y„ + y 3 ). 
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Now set F ( y ) = ln(-\/l + y 2 / y ) so that F '{ y ) = — (y + y 3 ) 1 - Then on (0,oo), F is 
positive, decreasing and convex, and we have 

(56) F ( y n+1 ) - F ( y n ) = F '^ n )( y n+1 - y n ) > > - CAt , 

yn + yi 

whence 

(57) F ( y n ) > F ( y 0 ) - CnAt = F { CM 0 ) - CnAt . 

Choosing any T* > 0 so that C* := F ( CM 0 ) — CT * > 0, we infer that as long as 
nAt < T* we have y n < U _1 (C*), and this together with lt5ll) yields the stability 
estimate ED- □ 

4. Unconditional stability for C 1 / C ° finite element methods without 

INF-SUP CONDITIONS 

The simplicity of the stability proof for the time-discrete scheme above allows us 
to easily establish the unconditional stability of corresponding fully discrete finite- 
element methods that use C 1 elements for the velocity field and C ° elements for 
pressure. 

We suppose that for some sequence of positive values of h approaching zero, 
Xh C H 2 fl 77 q (SI, K^) is a finite-dimensional space containing the approximate 
velocity field, and suppose Yh C i7 1 (S7)/R is a finite-dimensional space containing 
approximate pressures. We discretize ED in a straightforward way, implicitly only 
in the viscosity term and explicitly in the pressure and nonlinear terms. The re¬ 
sulting scheme was also derived in m and is equivalent to a space discretization 
of the scheme in m ED- Given the approximate velocity v!^ at the n-th time 
step, we determine p]^ £ Yh and u ^ +1 £ Xh by requiring 

(58) (VpJJ + pVV • Uh - vAu % + u ^- Vu ^ - f n , V0h) = 0 ^ h eY h , 

Vv n+1 — Vv 71 

(—^- - yv h ) + ( uAi % +1 , Av h ) = (Vp n h + - r, Av h ) 

(59) W h £ X h . 

We are to show the scheme above is unconditionally stable. First, we take 
< t>h = Ph in ED- Due to the fact that ( V(A — VV-)«J(, Vpjj) = 0, we directly 
deduce from the Cauchy-Schwarz inequality that 

(60) IIVKII < I|pVpsK)|| + K-V< - r II 

where 

( 61 ) VpsK) = (J-P)(A-VV-K 

is the Stokes pressure associated with u %. (Note Vp s ( uh ) need not lie in the space 
Yh). Now, taking Vh = Uh +1 in (1591) and arguing just as in ED: we obtain an exact 
analog of ED , namely 

^(||V< +1 || 2 - ||V<|| 2 ) + (u - £l )||A<+ 1 || 2 

(62) < ^ (ll/l 2 + K • V<|| 2 ) + H|Vp s «)|| 2 . 

Proceeding now exactly as in section 3 leads to the following result. 
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Theorem 3. Let LI be a bounded domain in M. N (N = 2 or 3) with C 3 boundary, 
and suppose Xh C H 2 (~l Hq(LI, R w ), Yh C 77 1 (fl)/R. Assume f £ L 2 ( 0, T ; L 2 (Ll, R^)) 
for some given T > 0 and vf^ £ Xh. Consider the finite-element scheme 
with iff]) - Then there exist positive constants T* and C 4 , such that whenever 
nAt <T*, we have 

n 

(63) sup \\Vu k h \\ 2 + Y^ IIAzi^|| 2 Ai < C 4 . 

0 </c<n , ,, 

— — Ac—0 

The constants T* and C 4 depend only upon LI, v and 

M 0h := ||Vu°J 2 + uAt\\Au ° h \\ 2 + f \\f\\ 2 . 

Jo 

In ClH further arguments are given that establish the convergence of these 
approximation schemes up to the maximal time of existence for the strong solution 
of the unconstrained Navier-Stokes formulation ®-®. 


5. Non-homogeneous side conditions 


Consider the Navier-Stokes equations with non-homogeneous boundary condi¬ 
tions and divergence constraint: 


(64) 

(65) 

( 66 ) 

(67) 


d t u + u ■ Vu + Vp = vAu + f 
V • u = h 


u = g 
u = u in 


(t > 0, x £ LI), 
(t > 0, x £ LI), 
(t > 0 , x £ r), 
(t = 0, x £ LI). 


What we have done before can be viewed as replacing the divergence constraint 
(l?)51) by decomposing the pressure via the formulae in © and © in such a way 
that the divergence constraint is enforced automatically. It turns out that in the 
non-homogeneous case a very similar procedure works. One can simply use the 
Helmholtz decomposition to identify Euler and Stokes pressure terms exactly as 
before via the formulae © and 0 , but in addition another term is needed in the 
total pressure to deal with the inhomogeneities. Equation © is replaced by 

(68) d t u + V(u-Vu — f — vA u) + Xp g h = uV(V • u). 

The equation that determines the inhomogeneous pressure p g h can be found by 
dotting with V(j> for <j> £ H 1 ^), formally integrating by parts and plugging in the 
side conditions: We require 

(69) (VPgh, V0) = -( d t (n • g), fi) r + (d t h, <j>) + (vXh, Vfi) 
for all <J) £ i7 1 (H). With this definition, we see from (I 68 I 1 that 

(70) (d t u,X4>) - { d t {n ■ g),4>) T + {dtfiA) = (iA7(V • u - h), S/fi) 
for every cf> £ i7 1 (H). This will mean w := V • u — h is a weak solution of 

(71) dtW = vAw in LI, n-Vw = 0 on T, 

with initial condition w = V • iti„ — Ai| t . So the divergence constraint will be 
enforced through exponential diffusive decay as before (see (EHI) below). 

The total pressure in now has the representation 

(72) p = p B + up s +Pgh, 
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where the Euler pressure p E and the Stokes pressure ps are determined exactly by 
© and © as before, and p g h is determined up to a constant by the forcing functions 
g and h through the weak-form pressure Poisson equation (©3) . Our unconstrained 
formulation of then takes the form 

(73) d t u + i-Vu + Vp E + pVp s + S7p g h = vAu + f (t > 0, x e 12), 

(74) u = g (t > 0, x e T), 

(75) u = uin (t = 0, ie!l). 

We shall state an existence and uniqueness result for strong solutions of the 
unconstrained formulation m ®. We refer to CEEI for the proof, which is based 
on using a classical trace theorem of Lions and Magenes to reduce the problem to 
one with homogeneous boundary and initial conditions. Then one uses the stability 
of the time-differencing scheme together with a standard weak compactness method 
to get the existence. 

Let 17 be a bounded, connected domain in R N (N = 2 or 3) with boundary L of 
class C 3 . We assume 

(76) u in e H\Ll,R N ), 

(77) f eL 2 (0,T-L 2 (Q,R N )), 

g e fi 3 / 4 (0,T; L 2 (r, R^)) n L 2 (0, T; H 3/2 (T, R^)) 

(78) n{g| d t (n-g) e L 2 (0,T; iT 1 / 2 ^))}, 

(79) he H h :=L 2 (0,T-,H 1 (Ll))nH 1 (0,T-,(H 1 )'(Ll)). 

Here (if 1 )' is the space dual to if 1 . We also make the compatibility assumptions 

(80) g = zi; n when t = 0 , x € T, 

(81) (d t (n-g),i) r = (d t h,l) Q . 

We define 

(82) V := L 2 (0,T;ff 2 (f7,K JV ))nJf 1 (0,T;L 2 (0,]R JV )), 
and note we have the embeddings (Hi p. 288], [El p. 176]) 

(83) V^C([ 0 ,T],H\n,R N )), H h ^C([ 0 ,T],L 2 (Q)). 

Theorem 4. Let LI be a bounded, connected domain in (N = 2 or 3) and 
assume \7b\ )- 1 Wit- Then there exists T* > 0 so that a unique strong solution of 


assu me j7b] ) 
1 731) - J75| ) e x\ 


exists on [0,T*], with 

u eL 2 ( 0, T*; H 2 (Ll, R N )) n H\ 0, T*; L 2 (12, R N )), 

V =v Vs + Pe + P g h € L 2 (0, T*; H 1 (LI)/R), 

where ps and ps are defined in © and W’ p g h e L 2 (H 1 ( 12)/R) satisfies 16 !A) . 
Moreover, ue C([0,T*],H 1 (LI,R N )) and 

v-u-he l 2 (o,t*-,h\li)) n h 1 { o,t* ; (h 1 )\li)) 

is a smooth solution of the heat equation for t > 0 with no-flux boundary conditions. 
The map fn || V • u — h || 2 is smooth for t > 0 and we have the dissipation identity 


“II V-Z-h \\ 2 + v\WV-u-. 


= 0 . 


(84) 
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If we further assume h G L 2 (0,T; H 2 (fl )) n H 1 (0, T ; L 2 {{ 1)) and V • zti n G 
then 


V • m G L 2 (0, T*;H 2 (n)) n ^(0, T*; L 2 (U)). 
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Appendix A. Uniform bounds on the Neumann-to-Dirichlet map 

Here we provide a more detailed sketch of the proof of inequality m in LemmaQ 
based on summarizing the relevant arguments from |lM 1 vol. I] and taking into 
account the finite regularity of the boundary. 

Without loss of generality we may assume fl is connected, so the same holds for 
Q s for small s (see below). For solutions of the Neumann problem 

(85) Ap = 0 in O,, ft-\7p = g on r s , 

we seek to bound the tangential gradient of p by g in L 2 (P S ), uniformly for small 
s. This is accomplished by interpolating between maps 

pr 1/2 (r s ) H\n s ) ^ ff 1/2 (r s ), 
iL 1/2 (r s ) ^ H 2 (n s ) % h 3/2 (t s ). 

Here the maps Tf and T 2 are the trace maps. The map A\ gives the zero-average 
weak solution of the Neumann problem by applying the Lax-Milgram lemma 
to the weak form of (1551) : 

(86) [ Vp-V<? = I gq for all q G H 1 ^^) with [ q = 0. 

Jq.% Jr 3 J n c 3 

For our given solution of <E3 we have p — p s = Aig where p s = 1 f Qc p. Since 


(87) 



lls , llff- 1 /2(r s )||T , i'7||ffi/2(r f ,) < C , o||fi'lli?- 1 /2(r 3 )lkllH 1 (fij) 


by a trace theorem IlMI vol. I, p. 41], the bound 

( 88 ) Ib-PslluRflj) < Clllffllff-Va^,) 

follows by taking q = p — p s and using Poincare’s inequality. 

The map A 2 is the restriction of A\ and is bounded by elliptic regularity theory 
using only the C 2 regularity of the boundary. Using that the trace maps T) and X 2 
are bounded, we obtain bounds 


(89) lb —-Psllffi/2(r„) < ^ 2 11S'11 1 / 2 (r s ) 5 lb _ P«llff 3 / 2 (r s ) < ^31 |^IItti/ 2(r s )• 

By an abstract interpolation theorem IlMI vol. I, p. 27], the map TiAi restricts to 
a bounded map between interpolation spaces, 


(90) TA : [H- V 2 ^), ff 1/2 (r s )] 1/2 - [tf 1/2 ( T„), H 3 / 2 (T s )} 1/2 
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with ||Tj4|| < C 4 where C 4 depends only on C 2 and C 3 . By |LMI vol. I, p. 36] we 
have the isomorphisms 

(91) L 2 (T s ) £* [HV 2 {T a ),H 3 ' 2 {T a )] 1/a , H 1 ^,) S [H~ 1 / 2 (T S ) 1 H 1 ^ 2 (T s )\ 1 / 2 

with equivalent norms, and so the map g 1 —> TAg = (p — p s )|r s is bounded from 
L 2 ( r s ) to i? 1 (r s ). This leads to the estimate 

II^Plli 2 (r s ) = ll^(P — Ps)lll 2 (r s ) 

(92) < ||p - Ps||/ri(T.) + II” '^ 7 PllL 2 (r a ) < Ci\\n ■ Vp||f,2(r s ), 
which corresponds to <GS>- 

We now argue that these bounds are uniform in s for 0 < s < so small. For 
this purpose one should consider how the fractional-order Sobolev spaces H r (T s ) 
and the norms in those spaces are defined, and how one obtains the bounds in 
trace theorems, in Poincare’s inequality, in the elliptic regularity estimate, and in 
comparing equivalent norms. All this can be done (except for Poincare’s inequality, 
which we discuss below) using a fixed partition of unity on f l and a fixed family of 
maps that locally flatten the boundaries T s simultaneously. 

To describe how this can work, let <P(x) = ±dist(x, P) be the signed distance 
function on R w , positive in Q and negative outside. Then for some so > 0, $ is C 2 
for 0 < |$(x)| < s 0 . Let Xq £ T be arbitrary. Translate and rotate coordinates as 
necessary so Xq = 0 and the normal n{x 0) = —e/v . For x = (x, xjv) £ x R, 

let 

(93) %/j(x) = (x,®(x,x N )). 

Evidently the map if) locally flattens all F s simultaneously — we have tl>(x ) = (x, s ) 
for all x £ r s . By the proof of the inverse function theorem, ip is a C 2 diffeomor- 
phism of i?(xo,ro) onto its image, for some tq £ ( 0 , so] fixed and independent of 
Xq. 

Let rfe = 2 -fc ro, k = 1,2, — The compact set Sl r2 is covered by the union of 
balls B(xo,ri) for Xo G F. Hence there is a finite subcover by balls B{x^\r\), 
j = 1,... ,m. We denote by i/jj the associated diffeomorphism on Bj = B{x^\r 0 ); 
we have efy • ipj(x) = $(x) for \x — x^\ < tq. We fix a partition of unity on VL r2 
subordinate to this finite cover by Bj. I.e., we fix ay £ Cg°(Bj), j = 1,..., m, such 
that X)JLi a j( x ) = 1 f° r x e &r 2 - 

The norms in the fractional-order Sobolev spaces H r (T s ) can then be defined 
as follows for |r| < 2. Given s £ [0,r 2 ], the norm of any u £ C 2 (r s ) in H r (T s ) is 
given by the norm of the C 2 function x i-> {(aju) o in Ff r (R Ar ~ 1 ) m . 

With this fixed system of cutoffs and flattening maps, one can check that the 
trace theorem for T\ and T 2 , the elliptic regularity theory for A 2 , and the norm 
equivalences that must be invoked to obtain the estimates in all involve bounds 
that are valid uniformly for s £ (0, 7 * 1 ]. It remains to check Poincare’s inequality. 

Uniform bounds in Poincare’s inequality. In the procedure that we have out¬ 
lined, to obtain a uniform bound in 1551) we need that the constant in Poincare’s 
inequality 

j \p-p,?< c J A = j4/» 


( 94 ) 
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can be taken independent of s for small s. Showing this is somewhat related to 
showing that f2° is connected for small s, which we can do as follows: Since 12 is 
connected, any two points in 12° can be connected by a path in f2. On this path 
we can replace points in f l s by the projection along normals onto IV, where s' — s 
is positive and so small that the new path is contained in 12 ° and still connects the 
same two points. 

To prove the uniform Poincare inequality (PH) , we reduce to a fixed s by integrat¬ 
ing along normals as follows. Recall that the distance function = dist(a:,r) is 
C 2 for <!>(a;) < so- Fix £ = 5 so, and suppose 0 < s < e. Let q be a smooth function 
on 12°. Then if s < d > (x) < e we integrate d t q(x — tn(x)) from t = 0 to e and apply 
the Cauchy-Schwarz inequality to get 

(95) \ q ( x )\ 2 < 2 \ q(x — en ( x ))\ 2 + 2 e f \Vq(x — tn(x))\ 2 dt. 

Jo 

Integrating this in the domain where s < < e and changing variables on the 

right hand side we get 

(96) f \q(x)\ 2 dx < C ( f W(y)\ 2 dy + [ \X7q(y)\ 2 dy 

Js<®(x)<e yj e<^(y)< 2e J s<®(y)<2c 

with C independent of q and s. Now, given a smooth function p on 12°, use this 
inequality for q = p — p s , together with the standard Poincare inequality for the 
domain 12 °, and the L 2 -optimality of the average in 12 °, to conclude that 

(97) / \p-p s \ 2 <[ \p-p e \ 2 <c([ \p-p e \ 2 +[ \Vp\ 2 )<c[ |Vp | 2 

Jqc \^e J ^ C s 

for some C independent of p and s. 
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